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Abstract 

In this paper, we put forward a way to study the nucleon's thermodynamic properties such as 
its temperature, entropy and so on, without inputting any free parameters by human hand, even 
the nucleon's mass and radius. First we use the Lagrangian density of the quark gluon coupling 
^ \ fields to deduce the Dirac Equation of the quarks confined in the gluon fields. By boundary 
conditions we solve the wave functions and energy eigenvalues of the quarks, and thus get energy- 
momentum tensor, nucleon mass, and density of states. Then we utilize a hybrid grand canonical 
ensemble, to generate the temperature and chemical potentials of quarks, antiquarks of three 
flovars by the four conservation laws of the energy and the valence quark numbers, after which, 
all other thermodynamic properties are known. The only seemed free paremeter, the nucleon 
radius is finally determined by the grand potential minimal principle. 
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In our Quark Gluon Coupling Model, the nucleons are described as confining the 
quarks (q) and gluons (g) inside them, which are the quanta of the fields. The quarks 
interact through the exchange of gluons, and the gluon couples to the conserved quark- 
current through gipjuipA 11 , where g is the coupling constant. The Lagrangian density for 
this model is 

C = ? [V K - gA,) - m q ] ip - -F, V F^ (1) 

where = d^A u - d u A^. 

Lagrange's equations yield the the Dirac equation with the vector field[l,2]: 

(id,* - gA,) -m q ]ip = 0, (2) 

and the equation of conserved quark current : 

d,F^ = (3) 

The energy-momentum tensor is: 

= ^F Xa F Xa + v^Yd v ip + d v A x F x ^. (4) 
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Lagrange's equations ensure that this tensor is conserved and satisfies d^T^ = 0. It 
follows that the energy-momentum P u defined by 



(fxT 



0u 



(5) 



which is a constant of motion. 

We observe that at high quark density, the vector field operators can be replaced 
by their expectation values, which then serve as classical, condensed fields in which the 
quarks move, 



(6) 



For a static, uniform system, the quantities A are constants independent of x^. In the 
Mean Field Theory, the Lagrangian density is 

C = ^[^- 5 7%-mJV (7) 
Hence, the Dirac equation is linear, 

Vfy - 91° A - m q ] 4> = 0. (8) 
We seek normal-mode solutions of the form ^(x^) = t/j(r)e~ lEt . This leads to 

[-ia.V + gA + [3m q }^(r) = E^{f) (9) 

Consider the case k = — 1, which is the S\/2 level. The normalized quark wave function 
for a sphere of radius R is [3]: 



ip q (r,t) = Me 



-ie a t 



X 



V 



e q -gAg-m q 



jiU(tq-gA ) -m 2 r)) 



Xg 



4ir 



(10) 



where r = |f|, % g is the quark spinor and Af is the normalization constant [4]. e q is the 
quark energy eigenvalue. The density of quarks is readily calculated as 



where 



e<z - 9A0 ~ m q . 2 



e q - gA + m q 



X (V( e </ - 9Ao) 2 -m\r 



V , 
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r <R 
r> R. 



Thus, the density certainly does not vanish at r = R. Clearly, although the lower 
component is suppressed for small r, it does make a sizeable contribution near the surface 
of the nucleon. However, it is not the density, but i/ji/j should vanish at the boundary in 
the relativistic theory [5], 



e q - gA + m 



9 foW(e q -9A c 



r=R 



m q r) — - 



e n - gA - m q . 2 
Ji ( 



gA f - mlr) 



= 0. 
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This yields the energy states 



/ nir 2 

e q , n = 9Ao + Jrn? q + (—) , (11) 



where n — 1, 2, ... and the density of states is 

Then we consider the energy-momentum tensor for the following type, 

and is the familiar energy- momentum tensor for a free Dirac field = iif} / y' 1 d u ifj. 

The condition for overall energy and momentum conservation is that the divergence 
of the energy-momentum tensor should vanish, and this is certainly true for T^ u , as is 
easily proven from the free Dirac equation d^T^ = 0. 

However, the fact that the quarks move only inside the restricted region of space V 
leads to problems. Indeed, d^Oy = n^A s , where A s is a surface delta function A s = 
—n.d(8v) and n M = (0, f). In the static spherical case we find that A s is simply 8(r — R). 
Putting all these together we obtain d fM Ty U = iip^.nd v ipA S) and using the linear boudary 
condition, 

W = -I A s = -Pn»A s , 

Z s 

where P is the pressure exerted on the sphere's surface by the contained Dirac gas 

P=-\ n.£>W) 



To keep the energy-momentum conservation, we add an energy density term B6y to the 
Lagrangian density. Then (since T^ v involves Cg^ u ) the new energy-momentum tensor 
T% v has the form 

TjT = (T" u + Bg^)6 v . 
Therefore, the divergence of the energy-momentum tensor is 



d li Tr = (-P + B)n v A 8 , 

which will vanish if 

B = P = - 
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B = P = ~ n.d u {ipi>) . (13) 



Thus the nucleon's mass is given by 
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M 



= J dx 3 9 v T 00 + -nBR 3 . (14) 



Then we study the nucleon's thermodynamic properties, considering it as a many- 
body system at a given temperature (T). This is a hybrid system within radius R of 
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interacting massive quarks, antiquarks of 3 flavors (u, u, d, d, s, s) and massless gluons 
(g) , all of which can exchange the energy and the particle numbers[6]. In a statistical 
approach, the grand canonical partition function is given by 

Zq — Z vac Z u Zu Z d Z d Z s Z ? Z g , (15) 

where Z vac takes care of the temperature T — > limit, and we consider this as 

-T\nZ vac = \itBR\ (16) 

Z u , Z d , Z s , Z u , Z d and Z s refer to the partition function for the u,d,s quark and u, 
d, s antiquarks, while Z g refers to the gluonic part. Their chemical potentials obey that 
ji(u) = —fj,(u), fj,(d) = —/i(d), fj,(s) = —fJ-(s) and fi{g) = 0, among which there are only 3 
free parameters. 

The total nubmer, energy and grand partition function of each quark system are: 



N i = Vl/f e ( ^"^ )/T + l)=6 r l/(e^-^l T + l)p q (e)de, (17) 

• 1 V ' Jm q +gA 1 V ' 

E q = E e 7( e^~^' T + l) = 6 [°° l/( e^' T + l)p q (e)ede, (18) 
= E ln ( 1 + e-^-^ T ) =6 f°° In ( 1 + e^-^We)*. (19) 

j Jm q +gA 

And the total energy and grand partition function of the gluon system are E g = 
^L-X 4 V and lnZ g = ^-T 3 V. We obtain the four parameters p(u), fi(d), fi(s) and T by 
the following four conservation laws [7]: 

N u -N u = n(u), 
N d -N d = n(d), 
N s - N- s = 0, 

and 

E u + E u + E d + E d + E s + E- s + E g = M - ^nBR 3 . 

Therefore the grand potential of the hybrid system is given by 

n = -T In Z G (20) 
Thus we get the entropy of the nucleon as: 

s = -(dn/dT)\ Ku)Hd)Hu)y (21) 

At last, the only unsettled parameter, the nucleon's radius R, that we input into the 
model at the beginning could be determined by the grand potential minimal principle: 

(dn/dR)\ Ku)Md)MT =0 (22) 
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